Non-analytic curvature contributions to solvation free energies: influence of drying 
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We investigate the solvation of a hard spherical cavity, of radius R, immersed in a fluid for 
which the interparticle forces are short ranged. For thermodynamic states lying close to the liquid 
binodal, where the chemical potential deviation 5fi = (i — ^i co (T) is very small and positive, complete 
wetting by gas (drying) occurs and two regimes of interfacial behavior can be identified. These 
are characterized by the length scale R c — /(Ap<5/i), where 7^ is the planar gas- liquid surface 
tension and Ap is the difference in coexisting densities at temperature T. For R> R c , the interfacial 
free energy and the density profile of the fluid near the hard wall can be expanded in powers of 
the curvature R~ , in keeping with the analysis of Stillinger and Cotter, J. Chem. Phys. 55, 3449 
(1971). In the other regime, R < R c , the interfacial free energy and its derivatives acquire terms 
depending on \nR. Since R^ 1 can be made arbitrarily small this implies non-analytic behavior, 
as R^ 1 — > 0, of the work of formation of a hard spherical cavity and of the Gibbs adsorption and 
the fluid density at contact with the wall. Our analysis, which is based on an effective interfacial 
Hamiltonian combined with exact statistical mechanical sum rules, is confirmed fully by the results 
of microscopic density functional calculations for a square- well fluid. We discuss the repercussions of 
our results for solvation phenomena, emphasizing that non-analytic behavior equivalent to that we 
find for complete drying in solvophobic systems will also arise in the case of complete wetting, i.e. 
when liquid films are adsorbed on the surface of large (colloidal) particles or at curved substrates. 
We re-assess various results in the important but neglected Stillinger-Cotter paper, where drying 
was not considered explicitly, in the light of our present analysis. 

PACS numbers: 68.08.Bc,05.70.Np 

I. INTRODUCTION 

If a solute is to dissolve in a solvent the latter must pay the price of forming cavities capable of accommodating the 
solute molecules. This process can usefully be imagined as being carried out in two steps, with the first step involving 
the insertion of a hard cavity in which all attractive solute-solvent interactions are turned off. Treating the basic 
model of hard cavity solvation is the aim of scaled particle theory (SPT)i One of the key issues in the chemistry and 
physics of solvation is understanding the significance of solute geometry and, in particular, the size-dependence of the 
free-energy of cavity insertion. If one naively equates cavity formation with bubble formation, then this issue reduces 
to understanding the curvature dependence of the surface tension of a bubble. Gibbs^ and later Tolman 3 . proposed 
analytic expansions of the interfacial free-energy in powers of the interfacial curvatures. This approach was adopted 
by the founders of SPT and remains today a commonly-used ansatz for treating inhomogeneous fluid phenomena. 

There is, however, one aspect of hard cavity-solvation that SPT is unable to describe even qualitatively. Namely, 
when the radius of curvature is large and the liquid solvent approaches liquid-gas coexistence, i.e. the chemical 
potential fi — > (T). In this special case, it is known from the theory of wetting/drying transitions that a thick 
film of gas must appear at the surface of the hard cavity. More specifically, the planar hard wall-liquid interface is 
dry (completely wet by gas) in the limit fi — ► /Lt+,(T) for all temperatures T where liquid-gas coexistence occurs^ 
It is straightforward to show that standard SPT is unable to accommodate the drying of a hard cavity in the large 
radius limit £ Recently, Chandler and co-workers have proposed that a full understanding of the hydrophobic effect 
(of much importance in biology as well as chemistry) requires the incorporation of drying phenomena as the size 
of hydrophobic solutes increased One of the key aspects of their interpretation is the realization that at standard 
temperature and pressure water is very close to saturation (bubble formation), i.e. the relevant thermodynamic field 
(chemical potential or bulk pressure) is very close to its value at bulk coexistence. Because of the ubiquitous presence 
of attractive power-law (dispersion) interactions between solute and solvent true drying films are not easy to create 
in the physical world. Even weakly attractive solute-solvent interactions will lead to partial rather than complete 
drying so that the thickness of the region of depleted density is greatly reduced. Nevertheless, it is argued that the 
incipient presence of drying should play an important role in hydrophobic phenomena at large length scales^ 
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The breakdown of SPT close to saturation implies more than a simple failure to predict the growth of drying 
films. Drying at a planar substrate is an example of interfacial critical phenomena and as such necessarily requires 
the existence of non-analytic contributions in, say, p — p co (T) to the interfacial free energy. For a curved substrate 
one might expect that this in turn implies non-analytic terms in the curvature expansions of surface tension and 
of surface order parameters (derivatives of the interfacial free energy with respect to thermodynamic fields). Such 
a scenario calls into question the whole basis of the Tolman and SPT appro ach . In fact, the formal status of the 
analytic curvature expansion has been controversial in the past, see e.g. Refs. I8l9i Unbeknown to those authors and 
apparently to most of the inhomogeneous fluids community, Stillinger and Cotter— had, in 1971, provided what, at 
first sight, appears to be a proof that for a general solvent the curvature expansion of the interfacial free energy 
is analytic in Rr 1 for spherical hard-body cavities of radius R. Notwithstanding the existence of this proof, the 
incorporation of wetting/drying phenomena at curved substrates has led authors to identify non-analytic In R factors 
in the curvature expansion of the interfacial free energy and order parameters, such as the adsorption and contact 
density, for solvents where the interparticle potential is finite-ranged^iii2ii£ The main purpose of our present paper 
is to elucidate the non-analytic behavior identified explicitly in Ref. ^3 and place it in the context of the proof of 
Stillinger and Cotter. Although the mathematical details and hence our discussion pertain directly to complete 
drying at hard spherical cavities, we emphasize that the non-analyticities and the cross-over behavior which we 
describe below are equally relevant to the physically important case of wetting films of liquid adsorbed on the surfaces 
of curved substrates and colloidal particles^ 

We focus, in the main, on model solvents where the interparticle potential is finite-ranged, say a square-well or a 
truncated Lennard- Jones fluid. The presence of dispersion interactions leads to different non-analyticities which will 
be discussed elsewhere^ Our starting point will be the phenomenological interface Hamiltonian appropriate to the 
wetting/drying of curved substrates, assuming finite-range interactions, 16 in which the binding potential is expressed 
as the excess Grand potential of fluid surrounding a hard spherical cavity of radius R. We consider a reservoir of 
liquid at chemical potential p > p co (T) whose pressure is denoted by p. If Vtot is the total system volume including 
the cavity then the volume accessible to the fluid is V acc — Vtot — 47ri? 3 /3. The wall-fluid boundary of infinite repulsion 
is located at r = R and this induces a film of gas to grow out to a radius R + £ where a fluctuating gas- liquid interface 
forms with surface tension j g i. The gas film is considered to have a pressure p gi the pressure of the metastable bulk 
gas at the same chemical potential p, so that its grand potential is —p g V g , where V g = 4ir[(R + £) 3 — R 3 ]/3 is the 
volume of the film. It follows that the excess grand potential can be expressed as 

Sl ex (£;T,p,R) ee n+ P v acc 

= {p-p g )V g +Ait[R 2 lwg + {R + i) 2 lgl +R 2 a{T)Ze* V {-£/!i}} , (1) 

where j wg is the surface tension of the wall-gas interface. The final term in Eq. |JT]|. involving the gas-phase bulk 
correlation length £ (the decay length of the gas tail of the planar gas- liquid interface) , is the leading-order mean- field 
interaction between the gas-liquid interface and the wall-gas interface. For a hard wall a(T) > for all temperatures 
T at which gas-liquid coexistence occurs. From the Gibbs-Duhem equation at fixed temperature T we can write 
(defining Ap = pi — p g , the difference of the coexisting densities) 

p-p g = ApSp , (2) 

where 5p = p — p co (T) denotes the difference between the chemical potential and its value at bulk saturation. The 
equilibrium thickness of the gas film £ eq (T, p, R) follows by minimizing the Grand potential as a function of I at fixed 
T, p, R (and hence at fixed pressures and surface tensions). To leading-order in £ eq / R this yields 

where superscript 00 denotes the planar limit at saturation. Noting that at worst £ eq is no bigger than order lni?, 
we can substitute (0 back into Eq. to identify the first two terms in a curvature expansion of the equilibrium 
interfacial free-energy: 

n™(R, p) = 4nR 2 { lwg (R, p) + lgl {R) + (£ + £ eq (R, p)) 

Hereinafter we suppress T, since all of our analysis will be isothermal, and in addition for any specific quantity 
we list only those variables that contribute at the order of the expansions discussed below. We also refrain from 
constantly reminding the reader of the presence of higher-order terms. To go beyond our analysis would first require 
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one to reassess the possible need to include additional terms and/or thermodynamic field variations of the quantities 
appearing in the interface binding potential Eq. For example, for our purposes we do not need to distinguish 
between ^ g i(R + £ eq ,p) and J g i(R)- 

The key conclusion of our work is already apparent from the introductory analysis above. Namely, the presence or 
absence of non-analytic In R contributions is controlled by a length scale 

27°? 

through the asymptotic approach to the planar limit at gas- liquid coexistence (where both R and R c are infinite), 

<S_J ln {#} (6) 

The length scale (JSJ) is identical to the length that controls capillary evaporation between planar hard walls, or, 
indeed, capillary condensation between planar walls that are completely wet' Chandler and co-workers^ have linked 
this length scale to the hydrophobic attraction between large solutes in water; for water at room temperature and 
pressure R c ~ 1.4^m, a surprisingly large length scale. 

In the following section we re-analyze the arguments of Stillinger and Cotter, that might otherwise lead the reader 
into rejecting outright the presence of non-analyticity (terms involving lni?) in the curvature expansion of interfacial 
free energies and order parameters" In Sees. III-V we detail the cross-over between a non-analytic and an analytic 
approach to the planar limit, both for the interfacial free energy and for interfacial order parameters. Sec. VI presents 
Density Functional Theory (DFT) results for a square-well fluid, that illustrate the existence of this cross-over. We 
conclude in Sec. VII with a discussion of the physical significance of our analysis and comments on some of the 
additional interesting material to be found in the paper by Stillinger and Cotter. 

II. STILLINGER AND COTTER SUM RULE ANALYSIS 

There are three exact results that can be applied to a spherical hard cavity immersed in an arbitrary solvent: 

4TrR 2 k B T [p w (R, (i) - ^(oo, p)] , (7) 



dnr q (R,fi) 



OR 

dn e *(R,u) r°° 

<9ju Jr 
^clR ^ = -4:irR 2 k B T-^-[p w (R,p)-p w (oo,p)] . (9) 

The first two of these results are virial route and compressibility route sum rules, respectively. A brief review of their 
derivation and physical content is given in the appendices to Ref. 0. The quantity p w (R, p) denotes the limiting value 
of the density profile p(r) at the surface of the cavity (r = R) and p w (oo, p,) = p/ksT, where p is the pressure of 
the reservoir, is its zero curvature (planar) limit. Eq. (JHJ) is the Gibbs adsorption equation for this problem; pb(p) 
is the reservoir density. The third equation © is the Maxwell relation that follows directly from the previous two 
sum rules. It is equivalent to the derivative w.r.t. R of Eq. (4.13) of Stillinger and Cotter' although the latter 
was derived by a much less direct route. Note that Eq. @ is interesting in that it relates the derivative w.r.t. the 
field p of a local quantity p w to the derivative w.r.t. R of an integrated quantity, the Gibbs adsorption T. Stillinger 
and Cotter consider the consequences of the assumption of an analytic curvature expansion of the adsorption T and 
the contact density p w , in the context of the Maxwell relation ©. We re- visit their argument. Let us introduce the 
variable z = r — R and assume following Ref. an analytic expansion of the density profile, thus 

/ N , X P'( Z ) 1 P"( Z ) 1 P'"(Z) , , 

p(z; R) = p(z; TO ) + ^Li + _^ + _^i, (10) 

p w (R,p) = p w (R) = Moo) + — — + 2j ~^2~ + 3| — i (11) 

where a prime denotes a partial derivative with respect to curvature (-R -1 ) and all the coefficients belong to the planar 
limit (at a given chemical potential). When (|10|) is substituted into the definition of the adsorption T, in Eq. JHJ), 
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we see that r(i?) will then also possess an analytic curvature expansion provided the integrals over the expansion 
coefficients exist (are non-divergent). In particular, the only obvious requirement is for the existence of the planar 
limit integral 

/•OO 

/ dzz 2 [p(z;oo) - Pb(p)} ■ 
Jo 

Interestingly, in the presence of an r~ 6 , but not r~ 7 , dispersion contribution to the intermolecular potential of the 
solvent, this integral is ill-defined because when an r -6 energy is integrated over the semi-infinite volume occupied by 
solvent it generates a z~ 3 asymptotic decay in the planar density profile p(z; oo) away from the wall^ Thus, one does 
in fact anticipate the presence of a term or terms involving lni? in the adsorption when r~ 6 dispersion interactions 
are included. Indeed it has been shown explicitly that the free-energy then has a lni? contribution? 14 - 20 However, 
this non-analyticity is of higher-order than the corresponding result for the wet ting/ drying case with finite ranged 
potentials detailed below. 

For a strictly finite ranged model potential one might expect the integrals over the expansion coefficients to exists 
and therefore the adsorption to possess an analytic expansion in powers of the curvature. From the compressibility 
sum rule (JSJ this property is also transferred to the interfacial free energy (apart from unlikely terms independent of 
the chemical potential). Assuming T(R) is analytic in i? -1 the right-side of the Maxwell relation 10 cannot contain 
a term of order i? _1 , so that dp'" (Q) / dp in expansion (f 1 1 1> must be zero for arbitrary values of the chemical potential 
p. Stillinger and Cotter then argue that for the case of vanishing bulk density, pt(p) — > 0, p"'(0) is zero and therefore 
p"'(0) = for all p. If this condition is not met then one has a clear inconsistency with @. Note that since there is 
no reason to expect p"'(z) to be zero away from the wall, a failure to appreciate the special significance of the Maxwell 
relation ||2J would lead one to expect a lni? term in the interfacial free energyi&2 In fact, the enforcement of this 
unexpected condition on the limiting value of the density profile at a hard wall cavity has recently been demonstrated 
numerically within density functional theory for a hard sphere solvent^ The corollary to the argument of Stillinger 
and Cotter is that non-analytic terms in the curvature expansion of the interfacial free energy demand analogous 
non-analytic terms to be present in the interfacial order parameters p w and Y. Although Stillinger and Cotter did 
not contemplate drying at a hard cavity it is clear that wetting/drying non-analyticity must behave consistently with 
their sum rule arguments. In the following sections we provide an explicit analysis, based on Eq. JIJ and the sum 
rules listed above, which ascertains leading order non-analytic behavior. 



III. SURFACE TENSION ROUTE 



Substituting Eq. © into Eq. we can write the interfacial free energy as 

n e e x g (R,p;) = 4wR 2 { 7wg (R, t x) + lg i{R) ( 1 + ^f) + £Ap5p 
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h U^ }. d2) 



ApSpR + 2% 
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The behavior of this quantity depends on whether R > R c or R < R c , where R c is defined by Eq. (JSJ. There are two 
qualitatively different regimes. Let us consider these two regimes separately, in each case discarding terms beyond 
the leading-order curvature dependence. 

For i? i? c , for which dp must be kept sufficiently small, we have to our desired order in the curvature expansion 

^(i?, M )=4^i? 2 | 7lug (i?, M )+ 7s T + ^ln|^|| . (13) 
If we now invoke sum rule Q) it follows that in this regime 

k B T[p w {R,n)-p w {™,p)] = -[ lwg {R,p)+ 1 %] — , (14) 



where we remind the reader that terms of higher-order than those specified are omitted. This is the result obtained 
by Evans et al. in Ref. 0, based on setting Sp = at the outset^ Here, we note that for any non-zero Sp the planar 
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limit must lie outside this regime; i.e. at R = R c there must be a cross-over to analytic behavior. 

For R 3> R c , an approach to the planar limit at non-zero 5p, the free-energy curvature expansion 1)12(1 changes 
character to yield an analytic form whose leading term is 



fi~(fl, ft) = AnR 2 lwg (oo, M ) + 7 3 T + &pSl* 



111 



[ApSpj 



When this result is inserted into sum rule (J7J one obtains 



k B T [p w (R, fx) - p w {oo, (j,)] = — 

H 



lwg(0O, + + £Ap5p 



In 



ApSp 



The next term is 0(R 2 ); one has an expansion analytic in R 1 . Eq. I(16|) can be re-expressed as 

2 

k B T[p w (R,p) - p w (oo,n)] = — 7 tt ,i(cx),/i) + 0(R~ 2 ) , 

K 



(15) 



(16) 



(17) 



where 7^(00, /x) is the surface tension of the planar wall-liquid interface at chemical potential p\ this quantity contains 
the SphiSp term associated with complete wetting/drying for short-ranged forces. Thus, provided R 3> R c , one can 
read off the planar surface tension from the leading term in the curvature expansion of the contact density (which is 
Stillinger and Cotter's route to their exact formula for the planar surface tension discussed in Sec. VII below) even 
when drying is present. 



IV. SURFACE ADSORPTION ROUTE 
When Eq. (|12|) is inserted into the compressibility sum rule iJSJ we obtain 

aR 



T(R,p)-T wg (R,p) = -47ri? 2 £Ap m . 



ApSfiR + 2% 



which can be re-written as 



T{R,p) - T wg (R,n) 
4ttR 2 



In 



1 1 

R + RZ 



b(T) 



(18) 



(19) 



where b(T) = ln(a/2^). One can expand the logarithmic term in Eq. (|19|) about lni? c and note that the radius of 
convergence (and hence the extent of the analytic regime) is R c /R < 1. One can also read off the limiting forms of 
the two regimes. For the non-analytic case at R R c , examined by Evans et al. in Ref. IT2I 



T(R, n) - T wg (R, n) -> -AirR^Ap In ■ 



and for the non-zero dp analytic case at R 3> R c 



T(R, p) - T wg (R, p) -> -47Tfl 2 ^Apln 



aR 
27^ 



Ap5p 



(20) 



(21) 



V. SURFACE MAXWELL RELATION ROUTE 



In this section we will check directly the consistency of our curvature expansions for the case of drying with the 
surface Maxwell relation For the analytic regime R 3> R c one can demonstrate consistency by simply inserting 
(116(1 into the right side of @ and comparing with the derivative of Eq. 1(21(1 . From both routes one obtains 



d(T - T. b 
dR 



inR^Ap In 



\Ap5p 



(22) 
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The non-analytic regime R <C R c is not so straightforward because one must first retain the leading order term 
varying with chemical potential (missing from the right-side of Eq. (11411 ). This requires us to keep the non-dominant 
term in the prefactor of the logarithmic contribution to the interfacial free-energy (|12|l : it is also helpful to expand 
the logarithm and note that at next to leading-order it cancels the sub-dominant £ApSp term: 



n%(R, p) - ^R 2 *yt wg {R, m) + 1T1 + £ 

Thus, keeping appropriate terms, at non-zero Sp we have 

2 



Ap5p 



27, 



gl 



ksT [p w (R, fi) - p w (oo,n)] = 



R 



R 



ApSp 



In • 



aR 



27, 



21 

R 



In ■ 



aR 

27f 



(23) 



(24) 



where despite appearances the prefactor of the logarithmic term is not missing a factor of two in front of 7^ . When 
Eq.()24p is inserted into the surface Maxwell equation © one immediately rederives the result for dT/dR obtained 
from (|20() . to leading order as required. 



VI. RESULTS FROM DENSITY FUNCTIONAL THEORY FOR A SQUARE- WELL FLUID 



The explicit results for surface thermodynamic functions that we have presented in earlier sections are based 
upon the ansatz Q for the excess grand potential of a fluid adsorbed on a hard spherical cavity. This ansatz, as in 
all interface Hamiltonian approaches to wetting/drying, relies upon the introduction of the length £, the thickness 
of the drying layer around the sphere. In a fully microscopic approach only the fluid-fluid pair potential should 
appear, once the wall-fluid (external) potential has been specified. The properties of the fluid are determined by the 
behavior of the average one-body density p{r), which depends, in turn, on the nature of the wall-fluid and fluid-fluid 
potentials, as well as on the thermodynamic state point (T, p) of the fluid reservoir. In order to test the predictions 
of the (coarse-grained) interface Hamiltonian approach we employ the same microscopic density functional theory 
(DFT) used in Ref. [T^ Specifically we consider an attractive square-well fluid adsorbed at the hard spherical cavity 
of radius R. 



The fluid-fluid potential is given by 



with e > and the wall-fluid potential is 



00 , r < a 

-e ,0 <r< 3ct/2 ' [Z0> 



The grand potential functional is taken to be 

Slv\p] = F id [p] + F ex [p] + [ drp(r)(V(r) - p) , (27) 



where Fid [p] is the Helmholtz free energy functional of the ideal gas and the excess (over ideal) part of the free energy 
functional is approximated by 

F ex [p] = F hs [p] + ljdrj dr'p(r)p(r')4> att (\r - r'|) , (28) 

with p(r) = p(r). Here Fh s [p] is the excess Helmholtz free energy functional of a hard-sphere fluid which we treat 
by means of Rosenfeld's^ successful fundamental measures theory. Attractive interactions are treated by means of 
a simple mean-field approximation^ we take the attractive part of the potential to be (p a tt(i~) = ~£ f° r r < 3cr/2 
and zero otherwise. The equilibrium density profile Peq(r) is obtained by solving numerically the Eulcr-Lagrange 
equation resulting from minimizing the functional (|27(l at fixed p, R, and T. As usual, fiy[p eg ] yields the estimate 
of the grand potential O and, hence, of the surface excess quantity defined in Eq. from which other surface 
thermodynamic quantities follow. As emphasized in Ref. Il2t this particular DFT approach has the advantages that 
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i) the coexisting densities pi and p g can be calculated precisely from the bulk free energy density arising from l|28[l . 

ii) the results from the DFT defined by (|27l) and (|28|) satisfy the Gibbs adsorption theorem (JSJ and the sum rule 
for the contact density p w (R, p) , (this was confirmed earlier in extensive numerical calculations for hard-spheres 
adsorbed at hard curved cavities^), and iii) the planar surface tension 7^ and the bulk correlation length £ of the 
wetting phase (in this case gas) can be obtained from independent calculations. Specifically, we compute 7^ from a 
separate DFT calculation for a planar free interface and we calculate £ by evaluating the leading-order pole of the 
fluid structure factor obtained by taking two functional derivatives of Eq. (|28[l . ~ Note that the quantity £ entering 
Eq. Q is also the true correlation length which determines the exponential decay of the bulk pair correlation function 
g(r) for r — » oc^i The functional (|28|l is, of course, mean-field in character in that it omits some of the effects of 
capillary wave fluctuations that occur in a wetting/drying film. 25 We shall return to this issue later. 

In our numerical work we chose to focus on the surface adsorption route to surface thermodynamic functions, i.e. 
we chose to test the validity of Eq. (|f 9|) within the context of the microscopic DFT approach. It should be clear from 
the previous discussion that if the form of ljI9|l is verified then all the relevant predictions from the coarse-grained 
approach must also hold within the microscopic treatment. Note that in Ref. ^2 we had already confirmed, via DFT, 
the validity of the predictions of the coarse-grained theory for dp = 0, i.e. for R c = 00. Here we are considering finite 
but large R c /cr. We must also focus on the situation where R/a ^> 1 since the mesoscopic arguments based on Q 
are reliable only when the equilibrium layer thickness l eq ja is large, i.e. when bp is small and the cavity radius R is 
large. 

In Fig.^we plot the difference in adsorption T(R, p) — T wg (R, p), divided by the surface area 4irR 2 , obtained from 
the DFT calculations, with T defined by the integral in Eq. JSJ), versus the dimcnsionless quantity x where 

x = -£a 2 Ap In / kBT ^ . ] (29) 

for fixed R c = 5000<r and ksT/e = 1. The predictions of Eq. (|19|l . namely that this plot should be a straight 
line with positive gradient unity, are satisfied accurately both for R <C R c and R ^> R c . We have confirmed that 
for R/a > 250 Eq. H19fl remains valid for a wide selection of values of R c /a. At smaller values of R deviations of 
order a/R become apparent to the eye. We conclude from the DFT results that provided R/a is sufficiently large 
(and Sp is sufficiently small) the adsorption is given accurately by Eq. ifHfl) and thus the excess grand potential is 
well-accounted for by the starting equation (0J. 

The results in Fig. ^ are appropriate to a typical colloidal particle in 'water' at I atmosphere. As mentioned in the 
Introduction, we take R c = 1.4/zm for water at room temperature and pressure. In Fig. [21 we chose to fix R c to be 
250(J, corresponding to a thermodynamic state much further from bulk coexistence, and varied the colloid radius R. 
The difference in adsorption varies linearly with x for large radii R but deviations from a straight line of slope unity 
can be discerned for smaller values of R. 

In Fig. |21 we plot the difference between the contact density at the spherical cavity and that a planar wall, i.e. 
[p w (R, p) — Pw{oo, p)], versus 1/R for a fixed value of 5p corresponding to R c = 250<r. For R> R c the data lie on a 
straight line whose slope is given by 27^(00, p)a 2 /ksT. Thus, the DFT results also confirm the validity of Eq. I|17|l : 
the planar wall-liquid surface tension 7,^(00, p) at a non-zero value of dp can be obtained from plots of the contact 
density versus 1/R provided one has sufficient data in the range R> R c . 

VII. DISCUSSION 

In this paper we have used an interface Hamiltonian appropriate to (mean-field) wetting/ drying to explore the 
curvature dependence of the interfacial free energy Q, e e ^(R,p). We have concentrated on the solvation of a hard 
spherical cavity of radius R, a problem tackled originally by scaled particle theory, since for this system the exact sum 
rules {ZUHlini) provide direct insight. In particular, the surface Maxwell relation © implies that terms non-analytic 
in l/R cannot be present in the curvature expansion of the interfacial free energy unless there exist corresponding 
non-analytic terms in the Gibbs adsorption T(R,p) and in the contact density p w {R 1 p). Our main conclusion is 
that in the neighborhood of complete drying one must distinguish two regimes of behavior: i) R <C R c , where 
non-analytic terms involving mi? (see Eqs. H13III4ll2*0T ) are present and ii) R 3> R c , where the curvature expansions 
involve only powers of 1/R. The length scale R c is given by Eq. JSJ; the closer the bulk liquid is to coexistence, the 
smaller is the chemical potential deviation dp and the larger is R c . In Ref. [12J we considered only the case R c = 00 
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(Sfi = 0), identifying the leading-order, non-analytic contributions and confirming their existence by means of DFT 
calculations. Here we elucidate the crossover between the two regimes and verify the predictions from the interface 
Hamiltonian by performing DFT calculations of the Gibbs adsorption for a wide range of R and R c - see Sec. VI. As 
the DFT yields density profiles and surface thermodynamic functions consistent with the sum rules I|7I9|) . it follows 
that the leading-order curvature contributions predicted by the interface Hamiltonian analysis are all consistent with 
our DFT results. 



Of course, both theoretical approaches omit effects of capillary-wave fluctuations. In three-dimensional systems 
with short-ranged wall-fluid and fluid- fluid forces the mean-field theory of complete wetting/drying is expected to 
be affected marginally by fluctuations; the upper critical dimension is d c = 3. The linear renormalization group is 
sufficient to handle these effects and amounts to a Gaussian smearing of the interfacial binding potential. 28 Since, 

to leading-order in l/R, incorporation of curvature merely replaces ApSp for a planar wall by (^Apdp + 2j^f / R^j 

we conjecture^ by analogy with results for the planar problem, that our present mean-field results for the leading 
non-analytic contribution to (i?, fi), the third term on the right of Eq. I|13|) . should be unaltered when fluctuations 
are included, apart from replacing the bulk correlation length £ by (1 +oj/2)£, for ui < 2m& Here u> = UbT / (Air^(; 2 ) 
is the usual parameter measuring the strength of capillary- wave fluctuations in d = 3. Note that in the case of 
power-law potentials, arising from dispersion forces, the upper critical dimension d c < 3 and fluctuation effects are 
not expected to alter the results of the corresponding mean-field treatment. 

We turn now to the physical relevance of our results. As mentioned in the Introduction, there are potential 
implications for understanding aspects of hydrophobicity at large length scales, specifically for big solvophobic solute 
particles and for planar substrates^ In order to appreciate some of these we recall from the definition in Eq. Q that 
the excess chemical potential for inserting a single hard cavity into the solvent at fixed (/i, T), equivalent to the work 
required to create an empty cavity of radius R, is given by 

l*%{R,lx)=p^KR i + nT q (R,ii) , (30) 

where p is the pressure of the reservoir. Thus, a theory for the interfacial free energy constitutes a theory for the 
excess chemical potential associated with the insertion of a hard-sphere into the solvent, precisely the quantity 
that SPT attempts to calculate. The present analysis shows that striking logarithmic contributions can occur in 
fit x s for large R and sufficiently small 5fi. It is clear that no simple extension of SPT can hope to accommodate 
such subtle contributions. Indeed drying per se is not incorporated into standard SPT treatments.^ Moreover, it 
is evident that theoretical approaches to solvation which start from an approximate description of the Helmholtz 
free energy of the homogeneous mixture and obtain fj,f^ s by taking the derivative with respect to the solute den- 
sity in the limit of vanishing solute will not normally be able to capture the subtle physics associated with drying films. 



It is likely that the physical situation for which our theory is most directly relevant is not drying but one in which a 
big colloidal particle is wet by a liquid film in the approach to liquid-gas (fi — > fi~ Q (T)) or liquid-liquid coexistence^ 
The arguments that lead to Eq. Q are valid for complete wetting by either fluid phase, provided the wall-fluid 
potential has finite range or decays on a length scale shorter than the bulk correlation length of the wetting phase. In 
these circumstances the interfacial free energy per unit area will acquire a non-analytic R' 1 lni? contribution provided 
R <Si R c = 2j'*i/(Ap\8fi\). Of course the sum rule is modified when the wall-fluid potential V w f(r) is not purely 

hard 8 i 29 4 30 



4tt / drr l p(r) — - inR'p (31) 
oR J dr 

However, we expect the density profile in the neighborhood of the wall to acquire contributions equivalent to those 
appearing on the r.h.s. of Eq. (|14|) . More specifically, for a complete wetting situation, with R <C R c , the first term 
on the right in Eq. (|31|l should acquire both a contribution 87ri?7^ and a non-analytic 8n^^f In R contribution 
from the fluid- fluid interface located near R + £ eq . In real fluids dispersion forces are ubiquitous and these give 
rise to power-law divergences of the wetting film thickness at a planar wall. Analysis of the appropriate interface 
Hamiltonians and extensive DFT calculations for a Lennard-Jones liquid adsorbed at a hard spherical cavity show 
that wetting/drying leads to non-analytic terms in the curvature expansion of the interfacial free energy and of the 
surface order parameters that are power laws in R^ 1 (not logarithms), provided R <C R c . For example, for R c = oo, 
the r -6 decay of the fluid-fluid pair potential gives rise to the adsorption increasing as i? 7 / 3 , a term proportional to 
R~ 2 / 3 i n the interfacial free energy per unit area, and a term i?~ 5 / 3 in the contact density^ 
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Throughout our discussion we have referred, somewhat loosely, to the presence of non-analytic contributions to 
surface thermodynamic functions. It is important to understand precisely what physical repercussions such contribu- 
tions, which involve lni?, might have. We choose to fix Sfj, to be small, so that drying films develop, but non-zero so 
that i? c is very large but finite, and vary the curvature R . For large curvature, R <C i? c , the interfacial free energy 
per unit area has a term in i? -1 lni?, (see Eq. JTSJ)- On reducing the curvature crossover occurs and for R^> R c the 
same quantity has an expansion in powers of R , i.e. 



AttR 2 



lw i{^,^) + 0(R- 1 ) . (32) 



It follows that the adsorption and the contact density differ from their planar limiting values by terms 0(R~ 1 ) and 
higher powers of the curvature; no singularity develops in derivatives as i?~ x — > + . This means that one can also 
consider the situation of negative curvature, i? -1 < 0, corresponding to adsorption from the liquid confined inside 
the hard cavity— (One can imagine solvent particles 'ghosted' into the cavity from a reservoir at the same fixed 
(//, T) - a situation also considered in Sec. VI of Ref. |T(J.) By continuity one expects thick films of gas to develop 
on the hard wall. Their equilibrium thickness is still given by Eq. @, and for the same films are thicker than 
at the planar wall since i? _1 < 0. Thus one expects a continuous (analytic) dependence of surface thermodynamic 
functions on curvature in the neighborhood of i? -1 = 0. As the curvature becomes more negative (the radius of 
the spherical cavity becomes smaller), the film thickens further and eventually the denominator of the logarithm in 
Eq. J2I vanishes when the curvature approaches — i?" 1 , signaling the breakdown of the theory. However, in reality 
this scenario will be prevented by the prior occurrence, at some larger radius, of capillary evaporation, i.e. the cavity 
will empty leaving only the 'gas phase—. 



If we make dfi smaller the regime R 3> i? c , where the surface thermodynamic functions exhibit power-law depen- 
dence on the curvature, shrinks. Nevertheless, provided R c remains finite there should be a (narrow) regime where 
Eq. (|32fl and the corresponding result (|17|) for the contact density remain valid. (The latter is especially important 
as it forms the starting point for deriving the planar wall-liquid tension in the analysis of Stillinger and Cotter— see 
Eq. Ull)).) Of course, one is free to work at bulk coexistence, R c — oo, and then the regime where a power series expan- 
sion in the curvature exists is of vanishing extent and T ~ — R 2 In R; this is the situation considered earlier in Ref. 
In this case capillary evaporation occurs for an infinitesimal negative curvature; the phase transition intervenes 
immediately to prevent our taking the logarithm of a negative number! Another situation one might contemplate is 
that in which the radius R of the hard cavity is fixed at some (large) value and the chemical potential \x is varied. 
For sufficiently large S/u,, R^>> R c and one expects power-law curvature contributions, whereas on approaching coex- 
istence Sfi — > and there will necessarily be crossover to the regime R >C R c where logarithmic contributions will arise. 

There are practical reasons for investigating the realm of validity of power-series expansions in the curvature. 
Helfrich's seminal analysis of fluctuating membranes is based on the assumption that the free energy can be written 
as such an expansion which includes three terms in addition to the bulk contribution— Very recently Konig et 
al— have argued that the interfacial free energy per unit area, W^(R, /j,)/ (AttR 2 ), and other surface thermodynamic 
quantities can be written as a constant plus only two contributions, one linear in the mean curvature and the other 
linear in the Gaussian curvature of the convex surface bounding the fluid. They support their conjecture with 
numerical DFT results for a hard-sphere fluid bounded by a curved wall. It is clear from our analysis that the 
occurrence of wetting/drying negates the possibility that power-series expansions, with a finite or an infinite number 
of terms, could provide a valid description of the interfacial free energy for every state point of the fluid. For fluids 
with short-ranged wall-fluid and fluid-fluid forces, exhibiting gas-liquid coexistence, the existence of a large interfacial 
length scale £ eq 3> £ precludes such a simple description. In the case of power-law (dispersion) forces the situation is 
worse. As mentioned in Sec. II, for r~ 6 interactions has a lni? contribution for a non-wet spherical cavity— *2S 

We return now to discussion of the Stillinger-Cotter paper— and enquire whether their conclusions require revision 
in the light of our present analysis. Recall that these authors postulate the existence of a power-series expansion <|f (J I) 
for the density profile and on this basis they conclude that the work of cavity formation, and therefore the interfacial 
free-energy, is "free of contributions varying logarithmically with the radius, in the large size limit." As we argued 
in Sec. II, their argument cannot be valid for power-law forces. In the case of short-ranged forces there were sound 
reasons for accepting initially the validity of expansions (|1U|> and Ijlll) but we have ascertained that the occurrence of 
drying may lead to logarithmic contributions to the contact density: compare Eqs. ljTT|) and (fHEjl . Consistency with 
the sum rules then requires the presence of logarithmic contributions in the interfacial free energy, in contradiction 
to Stillinger and Cotter's conclusion. However, if we restrict consideration to thermodynamic states and radii for 
which R ^> R c then the power-series expansions are appropriate and Stillinger and Cotter's conclusion could be 
deemed valid. We re-iterate that these authors did not contemplate drying so the length scale i? c does not enter 
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their analysis; effectively they set R c = a or £. 

Another important result of Stillinger and Cotter is their formula for 7^(00, fi), the surface tension of the planar 
hard-wall liquid interface - see Eq. (3.22) of Ref. Il0|. The formula, which involves integrals over the pair correlation 
function of the inhomogcncous fluid, can with hindsight be identified with exact compressibility route expressions for 
the planar surface tension (discussed below). Stillinger and Cotter obtain their formula by deriving a result for the 
quantity p'(0) entering Eq. i|ll|) . They then identify p'(0), the coefficient of the R^ 1 term in the curvature expansion 
of the contact density, with 27^(00, p)/kBT. This is equivalent to employing Eq. (jPo|) or (|T7Jl . The latter are valid 
provided R R c . However, this is not a handicap as one can always take the planar limit by fixing Sp to be small 
but non-zero and allow R^ 1 — > + . 35 As emphasized earlier, the resultant formula for 7,^(00, p) does capture the 
gas-liquid tension plus the 5p\u5p contribution which characterizes complete drying in systems with short-ranged 
forces. 

As with much of their paper, Stillinger and Cotter's derivation of the planar hard wall-liquid surface tension was 
ahead of its time and appears to have been missed or ignored by later workers. The subsequent history of their result 
is interesting since the formula played a key role in unraveling some of the more subtle aspects of wetting/drying. 30,36 
Usually the first published compressibility route expression for the planar surface tension is attributed to Triezenberg 
and Zwanzig^i It turns out that the earlier formula of Stillinger and Cotter is equivalent to the functional inverse of 
the Triezenberg-Zwanzig (pair direct correlation function) formula, with the latter applied to the hard wall interface 
rather than the free liquid-gas interface; see for example Appendix A of Ref.l^ It was not until some years later that 
Schofield 38 proved the equivalence of the two compressibility route expressions with the much earlier virial expression 
of Kirkwood and Buff. Subsequently, without being aware of the work of Stillinger and Cotter, one of us rederived 
their approach for the surface tension by taking curvature derivatives of the Grand potential in order to isolate surface 
terms in the curvature expansionist A succinct general derivation, valid for arbitrary dimensionality and wall-fluid 
potentials, is given in Sec. IIIC of Ref. Specializing to three-dimensions and hard spherical cavities, the general 
formula reduces to 

l kB T# ( 
2 \dRj^ 

k B Tp 2 w dr 12 rl 2 [g w (r 12 -R)-l] , (33) 
Jo 

where, as previously, -f(R,p) = Q,^{R, p)/{AirR 2 ) is the surface excess grand potential per unit area and 
p w = p w (R,p) is the contact density of the fluid. gw{ri2', R) is the pair distribution function for pairs of particles 
1 and 2 positioned at the surface of the cavity. It is clear that the exact result Q33JI allows one, in principle, to 
investigate curvature corrections to the planar tension; only the latter was extracted by Stillinger and Cotter. In 
earlier papers it was assumed that the leading curvature correction would correspond to a Tolman term. Now we can 
see that in the presence of wetting/drying in the R R c regime, curvature corrections must generate the stronger 
R^ 1 \nR term in j(R, p) predicted by Eq. (|13fl . 

The compressibility route sum rule (|33|l implies pronounced fluctuation effects in the presence of drying. These 
are well-documented for the case of a planar wall 30,36 but have not been investigated in any detail for a curved wall. 
Consider Eq. I133[I in the limits where R — > 00 and 8p — > 0, where the gas-liquid interface will be located far from 
the hard wall, near R + l eq . Nevertheless, ^{R, p) remains dominated by contributions from this distant interface. 
The inescapable conclusion is that both the one-body contact density p w (R — > 00, p) and the pair distribution, 
9w{ti2]R — > 00) contain much information about the distant gas-liquid interface; capillary- wave like fluctuations of 
the latter manifest themselves at the wall. In principle, knowledge of the capillary-wave correlation contribution 
to g w (r\2', R)2^ would enable one to derive both the gas-liquid surface tension contribution and the leading-order, 
R~ 1 hiR curvature contribution to ^{R,p) predicted by the interface Hamiltonian approach. This is an ambitious 
program, as previous work on the planar interface will confirm^! 

There are other interesting topics discussed in the Stillinger-Cotter paper. These include the free energy for 
a droplet with a fixed number of solvent particles and some speculations on the behavior of the interfacial free 
energy near the bulk critical point - see Sec. VII of Ref. 0. Both topics warrant further attention; we cannot 
do justice to these in the present paper. In concluding, however, we do wish to return to what we believe is an 
important and intriguing issue in the statistical mechanics of inhomogeneous fluids. This concerns the vanishing of 
the coefficient p"'(0) of the R~ 3 term in the expansion of the contact density p w - see Eq. <|11|) . As we pointed out 
in Sec. II, Stillinger and Cotter argue, on the grounds of consistency with the sum rules, that p"'{0) must vanish 
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at all thermodynamic state points^ - with the caveat (see Sec. VII of their paper) that their curvature expansions 
(|10l II II) will themselves become inapplicable at or near the gas-liquid critical point where the bulk correlation length 
is diverging. Stillinger and Cotter do not write down an explicit expression for p"'(0) but they do write down the 
corresponding expression (their Eq. (5.13)) for p'^iO), the coefficient of Br 2 in the expansion of the density profile 
of a fluid adsorbed at a hard disk in two dimensions. Their sum rule arguments lead them to conclude that p 2 'd(0) 
must be identically zero. Commenting on this result Stillinger and Cotter write: "The identical vanishing of p'^iO) 
in two dimensions is hardly a transparent property of the complicated p'^iz) expression." This is certainly an 
understatement. If one glances at their Eq. (5.13) one finds integrals over combinations of one, two and three-body 
distribution functions of the inhomogeneous fluid evaluated in the planar wall limit and it is not at all obvious 
why P2d{®) should vanish. The corresponding expression for p"'(0) in the three dimensional fluid is even more 
complicated, involving additional four-body distributions. It is remarkable that these complicated expressions should 
vanish identically but this appears to be a necessary consequence of enforcing consistency with the Maxwell relation . 

What then is the status of Stillinger and Cotter's conclusions regarding the vanishing of p"'(0) and /^(O)? 
Their arguments require the existence of the curvature expansion 1)10(1. This is invalid for power-law intcrmolccular 
potentials and must fail near the critical point for any type of interaction. We have shown that the expansion can also 
fail, owing to the onset of drying, for all states sufficiently close to gas-liquid coexistence. Given the existence of these 
exceptions, one might be concerned that the expressions derived by Stillinger and Cotter for p"'(0) and P2 d (0) in terms 
of integrals over correlation functions might not necessarily exist for all thermodynamic state points. However, one 
could restrict attention to states where the bulk correlation length £ <C R and where wetting/drying does not occur or 
to states where R > R c if wetting/drying does occur. In fact, this latter restriction is implied in footnote 8 of a slightly 
later paper by Stillinger 40 - another paper remarkable for its allusions to future developments; in particular, Stillinger 
almost arrives at a modern statement of drying (and hence wetting) transitions (see his Fig. 3 and footnote 5). 
Under these restrictions the expansions IjlQI lll|) should exist. Do we then accept that the Stillinger-Cotter sum rule, 
their Eq. (5.13), is ill-defined outside these restrictions? This is feasible because the right-side of their complicated 
compressibility route sum rule contains integrals over the planar distribution functions alone, which could become 
very long-ranged in the presence of capillary- wave fluctuations (or bulk critical fluctuations or power-law interactions). 

However, provided dp > 0, it is difficult to see why the sum of integrals should vanish for one state and then not 
exist for another close by. The correlation length for complete wetting diverges only at dp — 0. What is more likely 
is that the sum of integrals is zero for all state points with 6p > but in the non-analytic R < R c regime, where the 
power series expansion is not valid, there could be an additional, non-vanishing contribution to the R~ 3 term in Hll|) 
that implies a logarithmic term in the free energy^i 

These matters are not of purely academic interest. There are many practical applications, see e.g. Ref. 0> where 
it would be beneficial to have a theory that expresses the interfacial free energy and density profiles as power series 
expansions in the curvature. As mentioned earlier, DFT studies of hard-sphere fluids 22 provide compelling numerical 
evidence for the vanishing of p"'{0). Those studies are based on the Rosenfeld functional 2 ^ which has its basis in 
scaled particle theory, i.e. it is constructed from a finite number of fundamental (geometric) measures. However, it 
should not be the particular form of the Rosenfeld functional that dictates p"'(0) = 0. Any DFT satisfying the sum 
rules (|7I9|I must necessarily yield p"'(0) = at a hard cavity, provided the curvature expansion exists (in particular, 
R > R c ). If follows that the DFT must yield distribution functions that are consistent with the Stillinger-Cotter sum 
rule. The behavior of p"'(0) defined by our DFT applied to the non- analytic regime R < R c requires further analysis. 
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In the case of a hard cylinder of radius R the term in question is known to be ^^/(fesTJ? 2 ), for R < R c , which is certainly 
non-vanishing - see Eq. (8) of Ref. Il2l This does not, of course, prove that the right-side of Stillinger and Cotter's sum 
rule Eq. (5.13) remains well-defined. Note that for the sphere a R~ 3 term in the expansion of p w (R) implies a contribution 
proportional to In R in Qlg . This would be weaker than the R In R contribution associated with complete drying - see present 
Eq. ED. 
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FIG. 1: The difference in adsorption V(R,/i) — T wg (R, p), divided by AnR 2 , for a square-well fluid in the liquid and gas phases 
versus x = -£a 2 Ap \n{k B Ta~ 3 / (2^ (i? -1 + R^ 1 ))} for fixed temperature k B T/e = 1 and fixed R c = 2^/{Ap8p) = 5000cr, a 
value typical for water. The results of the DFT calculation (symbols and line) are for sphere radii ranging between R = 5 • 10 5 er 
(bottom) and R — 295<r (top). The line has gradient unity consistent with the prediction of Eq. 11191 . For this temperature 
Apa 3 = 0.54542, j^a 2 /k B T = 0.19327 and the correlation length of the coexisting gas phase £ = 0.48962cr. 



FIG. 2: As in Fig.Qbut now for fixed R c = 250a. The sphere radii range between R = 5 • 10 5 a (bottom) and R = 295a (top). 
Small deviations from a straight line with gradient unity can be ascertained at small values of R (large x). 



FIG. 3: Difference between the contact density at the spherical cavity and that at a planar wall, p w (R,p,) — p m (oo,p), versus 
a/R for a fixed R c = 250er and ksT/e = 1. The solid line denotes the DFT results. For R > R c these lie on the straight 
(dotted) line whose slope is given by the planar tension 27^(00, p,)a 2 /ksT, as predicted by Eq. Ijl7^ . 
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